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ABSTRACT 



Context. Baryon Acoustic Oscillations (BAO) are a feature imprinted in the density field by acoustic waves travelling in the plasma 
of the early universe. Their fixed scale can be used as a standard ruler to study the geometry of the universe. 
^ f Aims. BAO have been previously detected using correlation functions and power spectra of the galaxy distribution. In this work, we 

present a new method for the detection of the real-space structures associated with this feature. These baryon acoustic structures are 
spherical shells with a relatively small density contrast, surrounding high density central regions. 

Methods. We design a specific wavelet adapted to the search for shells, and exploit the physics of the process by making use of two 
diff'erent mass tracers, introducing a specific statistic to detect the BAO features. We show the eff'ect of the BAO signal in this new 
statistic when applied to the A - Cold Dark Matter (ACDM) model, using an analytical approximation to the transfer function.We 
confirm the reliability and stability of our method by using cosmological A/^-body simulations from the MareNostrum Institut de 
Ciencies de I'Espai (MICE). 

Results. We apply our method to the detection of BAO in a galaxy sample drawn from the Sloan Digital Sky Survey (SDSS). We use 
the 'Main' catalogue to trace the shells, and the Luminous Red Galaxies (ERG) as tracers of the high density central regions. Using 
this new method, we detect, with a high significance, that the LRGs in our sample are preferentially located close to the centres of 
shell-like structures in the density field, with characteristics similar to those expected from BAOs. We show that stacking selected 
shells, we can find their characteristic density profile. 

Conclusions. We have delineated a new feature of the cosmic web, the BAO shells. As these are real spatial structures, the BAO 
phenomenon can be studied in detail by examining those shells. 

Key words, cosmology: large-scale structure of Universe - cosmology: distance scale - galaxies: cluster: general - methods: data 
analysis - methods: statistical 

1, Introduction Baryon acoustic oscillations (BAO) are therefore due to 

J-H I the propagation of these s ound waves in the baryon-photon 



plasma in the early universe (Peebles &Yu' 19701; iHu et al.ll 19971: 

Before recombination, the energy of photons is high enough to Eisenstein & Hu 1998; Bassett & Hlozek 201()^ Any primordial 

avoid the formation of neutral hydrogen atoms. This means that over-density in the early universe produces a spherical acous- 

baryons and photons are coupled through Compton scattering tic wave in the baryon-photon plasma, travelling outwards: the 

and electromagnetic interaction between protons and electrons, radiation pressure drags the baryons that are coupled to the 

forming a plasma. In this fluid two phenomena act in opposite di- photons, and compensates the gravity force that puUs aU mat- 

rections: gravitational forces tend to compress the plasma around ter towards the centre. Dark matter, however, is totally decou- 

high density regions, while radiation pressure tends to dilute pled from the photons, and therefore its density at the centre 

any such over-density. The combination of both in the presence continues growing. About 380,000 years after the Big Bang, 

of any initial inhomogeneity give rise to acoustic waves prop- temperature drops so that photons and baryons decouple, and 

agating in the baryon-photon plasma. This phenomenon ends the scale of the baryon shells freezes. After this time, both the 

abruptly at the epoch of recombination, when the temperature central over-density and the shell grow gravitationally, accret- 

drops enough to allow hydrogen atoms to form, and therefore ing both dark matter and baryons. The result at late times is a 

radiation decouples from the baryons. large over-density at the position of the original perturbation, 
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surrounded by a faint s pherical shell at a fixed co-moving scale 
(lEisenstein et al.ll2007h . 

The BAO scale is fixed by the sound horizon at decou- 
pling: it is the distance that the expanding acoustic shells can 
travel before decoupling. It has been accurately measured by the 
study of the an isotropics in the Cosm ic Microwave Background 
(CMB) to be (iKomatsu et alJ l2Q09i) r, = 153.3 ± 2.0 Mpc = 
110.4 ± lAh~^ Mpc (v^here v^e take h = 0.72, Freedma n et al.l 
1200 llfl . Therefore, this scale, once measured, could be used as 
a standard ruler to measure the Hubble expansion rate wit h red- 
shift H{z) and the angular diameter distance D a (z) ( Coorav 'etaD 
120011: iBlake & Glazebrookll2003l; ISeo & Eisensteinl l2005). 

The BAO should appear as a series of damping wiggles in 
the matter power spectrum, with the locations of the peaks and 
throughs in ^-space being a function of and other cosmolog- 
ical parameters (Eisenstein & Hu 1998). All the harmonics sum 
up to the same peak in the galaxy correlation function ^(r) at 
the scale r^, and therefore it could seem more appropriate to use 
this statistic for the detection of the BAO feature on the avail- 
able galaxy redshift surveys encompassing large volumes of the 
universe ( Sanchez et al. 2008). 

The first detection (claiming a ~ 3 cr level) was reported in 



the analysis of the correlation function (lEisenstein et al.l 



of the Sloan Digital Sky Survey (^SD SS) (lYorketal 



2005 ) 



2000h 



Luminous Red Galaxies (LRG) sample (Eisenstein et al. 2001), 
and later in the power spectrum CCole et al. 2005) of the 2-degre e 
Field Galaxy Redshift Survey r2dFGRS) dCoUess etal.ll200lh 
But ce rtainly this is a controversial topic. 'Cab re & Gaztanagal 
(l201lh are not finding such level of detection using a data set 
twice as large in volume and in number of galaxies. They do 
not claim this result to be in contradiction with the standard 
ACDM model, but to b e a consequence of insufficient data. 
One of the arguments in ICabre & G aztanagal (1201 ll) is the fact 
that mixing model selection with parameter determination can 
lead to some confusion in the interpretation of the results and 
their significance. Diflferent authors are using diff'erent criteria 
to asses s the significance of ther BAO detection. For example, 
when (Eise nstein et al.l l2005b affirm that the baryon signature 
was detected at 3.4 cr (or at 3.0 cr when including only data 
points between 60 and 180 /z"^ Mpc) they are comparing their 
results of the SDSS-LRG correlation function with the expected 
for the best-fit pure CDM model an d different BAO mo dels. 
The best BAO detection up to now (Perci val et al.l 12010 ) was 
obtaine d studying the com bined power spectrum of LRG and 
'Main' iStrauss et al.l (l2002l) samples of SDSS, together with the 
2dFGRS sample, and is at the ~ 3.6cr level. The authors ex- 
plicitly state that since this number is obtained comparing to an 
arbitrary smooth model, the significance cannot be directly com- 
pared with the one reported in (Eisenstein et al. 2005). This is a 
clear example of diff'erent authors using diff'erent ways to assess 
the signific ance of their results that in practice are not compara- 
ble. iHiitsil ([2OO6) calculated the redshift space power spectrum 
of the SDSS-LRG sample drawn from the Data Release 4. He 
concludes that BAO models are favored by 3.3(T over the corre- 
sponding models without any oscillatory behavior in the power 

spectrunL 

Perc ival et al ] (I2007h detected BAOs in the clustering of the 
combined 2dFGRS and SDSS main galaxy samples, and use 
their measurements to constrain cosmological models, in partic- 
ular a given combination of the ang ular diameter distance Da(z) 
and the Hubble parameter H(z). Cabre & Gaztanaga' ('2009a" b) 
studied the LRGs anisotropic redshift space correlation function 



h is the Hubble constant in units of 100 km s ^ Mpc 



^((T, n), where n is the line-of- sight or radial separation and cr 
is the transverse separation. Moreover, 'Gaztanaga et alj (l2009l) 
have shown ho w to constrai n H(z) using the correlations in the 
radial direction. iKazin et al.l (1201 Oal) found similar results for the 
correlation measurements and uncertainties, but manifest dis- 
agreement in the interpretation of the results regarding the de- 
tection of a line-of- sight baryonic acoustic feature. 

More recent studies (Martinez etakl 120091: 
Cabre & Gaztanagal l2009al: ISanchez et al.l 120091: iKazin et al.l 
2010bh lave confirmed this detection in the last Data Release 
(DR7, lAbazaiian et aP l2009l) of the SDSS-LRG, containing 
twice as many galaxies as the original sample, although the 
observed peak is in these cases wider than that observed in the 
original detection - an issue that needs further explanation. 
These measurements of the BAO scale at a low redshift, com- 
bined with other cosmological probes, have been used to put 
stringent const raints on the values of cosm ological para meters 
(Tegmark et al] 120061 iPercival et al.ll20 07l: Isinchez et aP 120091 : 
Percival et al. 2010; Reid et al. 2010; Ka zin et al.ll2010ah . 

While ^Bassett & Afshordil (12010) argue that low-level detec- 
tions may not be sufficient to rob ustly estimate the cosmologi- 
cal parameters, ICabre & Gaztanagal (12011 ) show instead that it 
is still possible -assuming a model- to locate the BAO position 
with data providing very low significant BAO detection. 

It is important, therefore, to find evidence of BAO in the 
galaxy distribution based on complementary methods. A step 
further is to search for real structures in the galaxy distribution 
that are responsible for the BAO feature in these second order 
statistics. The detection of these structures would be a confir- 
mation of the existence of the baryon acoustic phenomenon. 
Moreover, if we are able to localize these structures in config- 
uration space, this would allow us to study in more detail the 
properties of BAO. 

In this paper, we introduce a new method for the detection 
of BAO, which is closely tied to the underlying physics of the 
process, and apply it to a sample drawn from the SDSS cata- 
logue. This method (described in Section |2]) is based on ana- 
lyzing directly the 3D galaxy distribution using a very specific 
wavelet function (which we called 'BAOlet'), which is specially 
well suited to search for BAO features. The method makes use 
of two diff'erent tracers, one to map the overall density field (in- 
cluding the BAO shells), and the other to locate the position of 
the largest overdensities, which should correspond to centres of 
the shells. As we study directly the galaxy distribution in config- 
uration space, this method also allows us to identify regions of 
space where the BAO signal is stronger or fainter. We describe 
the expected signal in the ACDM model in Section [3l using both 
analytical prediction and a A/^-body simulation catalogue. We de- 
scribe the samples used in the case of SDSS in Section IH In 
Section [51 we show the results obtained in this case. We also 
make a test to assess the significance of these results, and ex- 
plore the implications of this analysis regarding the localization 
of BAO structures. Finally, we summarize our conclusions and 
discuss possibilities for future work in Section [6l 

2. The wavelet detection method 

The basis of the new BAO detection method is to focus on the 
positions of massive dark matter haloes, which correspond to the 
location of large initial perturbations, and to look for the pres- 
ence of structures resembling the acoustic shells around these. 
Once we locate the positions of the large over-densities, we need 
to study the density field to identify the structures correspond- 
ing to the acoustic shells around these centres. An appropriate 
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method for iden tification of structures in continuous fields is 
wavelet analysis (Martine z et al.lll993l:lstarck & M urtagh"2006l: 
Ijones 2009). Wavelet transfor ms are widely used in m any areas, 
especially in image analysis (' MallatI l2008l: IStarck et alJ 12010). 
They are specially suited for the analysis of data at diff'er- 
ent scales, and identification of characteristic patterns or struc- 
tures. Wavelets have been used in Cosmology for the analy- 
sis of the large-scale s tructure, and of th e CMB aniso t ropics 
(iMartinez et aD Il993l: iRauzv et al.1 1 19931: IVielva et al.1 I2Q04I: 
IStarck et al.ll2006l: ISaajl2009^. 

Standard wavelet functions like the Mexican hat are, how- 
ever, not suitable for the detection of shells. Instead, we need 
a family of wavelets whose shape matches the type of struc- 
tures we want to find in our data. Therefore, we use a specially 
designed wavelet (the 'BAOlet'), well adapted to the search of 
BAO features - shell-like structures around our selected centres. 
We design this new family of wavelet func tions as a t ransfor- 
mation of the wide-used B- spline wavelets (ISaai]l2009l) . These 
i/^R,s(^) functions are spherically symmetric, and their radial pro- 
files are defined as 



x10 



r-R\ 



(1) 



where R and s are the two parameters that define the scale and 
width of the BAOlet function, aR^s is the normalization constant 
defined so that 

Ur,s\\^ ^ J l^/^,.(x)pdx =U (2) 

and Bsix) is the box spline of the third degree, defined by 

B3(x) = -^{\x - 2p - 4|jc - Ip -h 6|jcp - 4|x -h Ip |x 2p j . 

The BAOlet function is shown in Fig. [T] It can be thought of 
as a spherical shell of radius R and width s, with zero ampli- 
tude at its centre and therefore adapted to the detection of spher- 
ical shells of a given radius. This specific choice is motivated 
by the fact that the integrated profile is the widely-used one- 
dimensional 'B -spline' wavelet function that has a null mean 
and compact support [-2, 2]. These properties directly translate 
onto the BAOlet that has also a null mean -a requirement for any 
wavelet function- if R > 2s, and takes non-zero values only for 
R-2s <\s\<R + 2s. 

We describe the density field using the density contrast ^(x), 
defined as 



Po 



where p(x) is the density field, and po is its mean. Then, given 
a density contrast map ^(x), properly normalized as in equa- 
tion Q, we can construct, for each point in the parameter space 
(R, s), a BAOlet coefficient map as the convolution of our density 
field with the corresponding wavelet: 



■x)dV. 



(3) 
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Fig. 1. The BAOlet function. Here we show a 2D plot (bottom) 
of the wavelet <A/?,^(x) used in the analysis, as defined by equa- 
tion ([T]). The top panel shows a ID slice along the dashed-dotted 
axis. The wavelet is plotted here for R = 105 Mpc, s = 
30 Mpc. The red dot marks the centre of the wavelet. This 
function has a null mean (provided that R > 2s), and compact 
support. It takes non-zero values only for R - 2s < \x\ < R -\- 2s. 



is a very delicate aspect of the BAO detection in the two-point 
correlation function. 

Due to the properties of the wavelet, the coefficient maps 
Wr^s(^) should have a null mean when averaged over all points in 
the volume considered. Equivalently, if we sampled these maps 
at N random points uniformly distributed in the volume (x^^^), 
the expected value of the average of the coefficients is zero. 



e[{WrAx^%} = 0. 



(4) 



The BAOlet acts as a matched filter, which is sensitive to data 
containing shells of different radius and different widths. Its 
property of zero mean is also of high importance since it makes 
the statistics derived from the BAOlet coefficients independent 
of the background level. Indeed, it is obvious that any constant 
added to the input data would not change the BAOlet coeffi- 
cients. In comparison, the estimation of such a baseline level B(R, s) = (Wr^s(^^c^)) 



This condition holds even in the presence of shell-like structures 
in the density field. Of course, for such structures the value of 
WrA^c) (Xc is the centre of the shell) is positive, and remains 
positive in nearby points. For an ideal S(r - R) density shell the 
radius of the region around the centre where the wavelet ampli- 
tude is positive, is s; the positive signal in this region is compen- 
sated by negative amplitudes around |x| = R. However, if we are 
able to identify the positions of N massive haloes in the same 
volume (Xc^), we can define a new statistic B(R, s) as the mean 
value of the coefficients WrA^) at these positions: 



(5) 
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If there are indeed shell-like structures around the selected den- 
sity maxima x^^ as expected for baryon acoustic structures, we 
should find positive values of B(R, s) with the maximum of B at 
the (R, s) values characterizing these shells. 

We can obtain further information from the wavelet coeffi- 
cients Wr^s(^), as we have information on the actual dependence 
of the signal picked up by the BAOlet function on the position. In 
particular, fixing a set of parameters of interest (Ri, Si), we could 
use the coeflftcients Wr.^sMc) to identify which of the selected 
massive haloes are giving the largest signal for these character- 
istics of the shells. In the context of BAO, the parameters Ri, Si 
can be chosen a priori using a theoretical model, or a posteri- 
ori using the parameters for which the function B(R, s) attains 
it maximum. In this way, we can localize in configuration space 
the structures responsible for the largest BAO signal in a given 
sample. 

For our calculation of B(R, s), we sample the (R, s) parameter 
space on a grid. For each point (R, s), we calculate the coefficient 
map Wr^s(^) as the convolution of the BAOlet with the density 
field (equation O. We perform the convolution in Fourier space 
using a Fast Fourier Transform (FFT) technique. To avoid prob- 
lems with the FFT, we zero-pad a large region around our density 
cube. To obtain B(R, s), we sample Wr^sOO at the position of the 
N selected centres, and calculate the average value (equation O. 

Therefore, to apply this method, we need a way to map the 
overall density field ^(x), but also to locate the position of mas- 
sive matter haloes x^^ We have to use two diff'erent populations 
of mass tracers, so that they play the appropriate role in the de- 
tection algorithm. The idea of using two diff'erent tracer sets, one 
for the small perturbations and another for the high peaks, in 
a cros s-correlation analysis was anticipated by Eisenstei n et al.l 
(l2Q07 h. We implement here a similar idea, but using a wavelet 
tool directly on the density field. As detailed below, we use 
galaxies from the 'Main' and LRG samples of SDSS in this case. 
However, this choice would depend on the kind of data available 
in each case. 



3- Prediction from ACDM 

In order to better understand our method, we show here which 
results we expect according to the ACDM model, and the ef- 
fect of BAO in our new statistic B(R, s). We use for this aim 
both the analytical a pproximation to the transfer function of 
lEisenstein & Hul (119981 and the results from the MareNostrum 
Institut de Ciencies de I'Espai (MICE) simulation (iFosalba et al.l 
l2QQ8h . 

In the first place, we use the ACDM transfer function, which 
allow us to study directly the eff'ect of the BAO. However, in this 
case, we must do a series of approximations in order to make a 
prediction for B(R, s). We want to predict which is the typical 
result for the wavelet coeflftcient Wr^s at the position of massive 
matter haloes x^, as a function of R, s. From equation ([3]), we see 
that this is equivalent to study the typical density profile around 
such haloes, ^(y - x^). The ACDM transfer function allows us to 
calculate this profile, provided we know which is the initial per- 
turbation corresponding to the selected haloes. We make here 
the simple approximation of considering that these initial pertur- 
bations are point-like and spherically symmetric, and can thus 
be simply described by a Dirac delta function in configuration 
space. This corresponds to a constant value in Fourier space. As 
the transfer function T(k) describes the relative evolution of the 
diff'erent Fourier modes, the present day radial density profile 



cor responding to such ini tial perturbation will we given simply 
by (lEisenstein et al.ll2QQ7h 

p(r) = Cf (r) (6) 

where T(r) is the Fourier transform of the transfer function T(k), 
and C is a normalization constant that depends on the details of 
the initial perturbation, and on the cosmic growth function Di(z). 
From equations ([3]) and ([5]), we see that the eff'ect of C will be 
just to change the normalization of our statistic B(R, s). 

Wq used the fitting form ulae to the transfer function T(k) 
from lEisenstein & Hu (1998), and obtained the expected Wr^s at 
a large overdensity using equations (O and (O. In order to high- 
light the particular signature of BAO, we also calculated Wr^s 
using the 'no wiggle' transfer function formula, in which the 
BAO have been edited out. We used here the values Qm = 0.25, 
Qa = 0.75, Qh = 0.044, and h = 0.7 for the cosmological pa- 
rameters, to all ow a direct comparison with the MICE simula- 
tion. Following lEisenstein & Hul(ll998l) . the sound horizon scale 
in this case is = 109.3 /z"^ Mpc. The results for both cases 
are shown in Fig. [2l In the plot, we mask the region R < 2s, 
as for these values of the parameters our BAOlet is not compen- 
sated (its mean is diff'erent from 0). Comparing both panels of the 
Figure, we see clearly which is the eff'ect of the presence of BAO 
in our statistic. In the case without BAO Wr^s is always negative, 
and it presents a smooth gradient across the (R, s) plane. This 
gradient is due to the overall shape of the radial profile (equa- 
tion dp. However, in the presence of BAO, Wr^s shows a promi- 
nent peak with positive values. This clearly shows the idea be- 
hind the B(R, s) statistic. The BAOlet i//r^s acts as matched filter 
with a shape adapted to detect BAO shells. Therefore the positive 
values in the coeflftcients Wr^s correspond to the cases in which 
the radial profile is matched by the BAOlet shape. The values at 
which Wr^s attains its absolute maximum, T^^ax = HO/z"^ Mpc 
and ^max = 22 Mpc, correspond thus to the characteristics of 
the shell that best matches the observed profile about the selected 
centres. 

In order to test the reliability of the method, and of this 
ACDM prediction, we calculated the B(R, s) for a halo cata- 
logue drawn from the MICE simulations. We use d the publicly 
avail able halo catalogue from the 'MICE3072' run (ICrocce et al.l 
12010 ). This particular run contains 2048^ particles in a box of 
side 3072 Mpc, therefore covering a volume of 29 Gpc^. 
The simulation was run with the GADGET-2 code (SpringeJ 
12005), assuming a ACDM model with the parameters mentioned 
above. Haloes in the simulation were selected using a friends-of- 
friends (FoF) algorithm. 

We used the resulting halo catalogue at z = 0, which contains 
a total of 2819031 haloes containing 143 or more particles. This 
corresponds to haloes with masses > 3.35 x 10^^ Mq. The 
halo number density is thus 9.72 x 10"^ h^Mpc'^. We used the 
full halo catalogue as a tracer of the overall density field. We then 
selected as centres for the calculation of B(R, s) in equation ([5]) 
only the haloes with a mass > 1.76 x 10^"^ Mq. We chose this 
mass threshold in order to select approximately the 10% most 
massive haloes in the simulation box. This choice is somewhat 
arbitrary, but serves for the purpose of testing the BAOlet method 
and illustrating the expected result. 

Fig. [3] shows the BAOlet result B(R, s) for these MICE sam- 
ple s, compared to the theo retical results obtained above from 
the Eise nstein & Hul (1 19981) transfer functions. We obtain a re- 
sult very similar to that of Fig. [21 as B(R, s) shows a clear 
peak, and attains its absolute maximum for T^^ax = 108 /z"^ Mpc, 
Smax = 28 /z"^ Mpc. This indicates that our BAOlet method can 
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Fig. 2. Values of the BAOlet coefficients Wr^s expected at the 
positions of large initial point-like perturbations, as a function 
of the BAOlet parameters (R, s). The bottom panel shows the 
result using a standard ACDM transfer function, while the top 
panel shows the result using a transfer function with the BAO 
wiggles edited out (Eisenstein & Hul l 19981) . The normalization 
is arbitrary. The contours are drawn at steps of 1000 for Wr^s < 
(dotted), Wr,s = (solid), and Wr,s > (dashed). The map 
attains a maximum atR = llOh~^ Mpc, s = 22h~^ Mpc. 



be applied to two sets of mass tracers, although the details of the 
tracers used here are very different from the ones we use later 
on the SDSS samples. This also confirms the expected eff'ect of 
the presence of BAO in the B(R, s) function: the presence of a 
large peak with positive values of B, located approximately at 
the values of R and s corresponding to the radius and width of 
the acoustic shells. The fact that we obtain here slightly diff'erent 
values for T^max and ^max than those predicted above may be due 
to non-linear evolution eff'ects, which slightly reduce the radius 
and increase the width of the she lls. A similar effect is present in 
the correlation function (see e.g. lCrocce & Scoccimarroll2008l) . 



10 I 1 1 1 1 1 1 1 r 




-40 -30 -20 -10 10 

B(R,s) 



Fig. 3. The BAOlet statistic B calculated for the MICE simula- 
tion sample described in the text as a function of the parameters 
(R, s) (bottom panel). The contours are drawn at steps of 5 for 
B < (dotted), B = (solid), and 5 > (dashed). This func- 
tion attains its maximum for R = 108 /z"^ Mpc, s = 2Sh~^ Mpc. 
The top two panels show cuts at the values s = 2^h~^ Mpc (top) 
and s = 22 Mpc (middle), marked with grey horizontal lines 
in the 2D panel. In each case, the solid blue line corresponds to 
the value obtained from MICE, the dashed red line corresponds 
to the theoretical expectation from the Eisenstein & Hu (1998) 
transfer function (bottom panel of Fig. O, and the dotted green 
line to the theoretical expectation using the 'no wiggle' trans- 
fer function (top panel of Fig. [2]). These theoretical predictions 
have been re-normalised to get the same value at the maximum 
in B(R, s). 



We also used this halo catalogue from MICE to make a qual- 
itative estimation of how diff'erent observational eff'ects would 
aff'ect the BAOlet result. In the first place, we studied the ef- 
fect of redshift- space distortions. To this end, we calculated the 
redshift- space positions of all haloes taking into account their 
peculiar velocities, as output by the simulation, and consider- 
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ing an observer located in one of the vertices of the simulation 
cube. The result for B(R, s) in this case is shown in the top panel 
of Fig. m where it is compared to the real-space result discussed 
above. As can be seen from the figure, although small diff'erences 
appear between the real- and redshift- space results, the main fea- 
tures of the B(R, s) prediction remain the same, with the position 
of the maximum changing only by ~ lh~^ Mpc. 

For the second case, we added the eff'ect of a decreasing ra- 
dial selection function across the sample. We model this selec- 
tion as an exponential decay function, such that the final num- 
ber of haloes used to trace the overall density field is ~ 20% 
of the total. In our calculations, we then weight each halo by 
the inverse of the mean density at its redshift, as we do later for 
the SDSS data. We do not apply any selection function to the 
centres. The results for B(R, s) obtained in this case (including 
also the redshift- space eff'ects) are shown in the bottom panel of 
Fig. m As above, these observational eff'ects do not change sig- 
nificantly the overall behaviour of B(R, s), or the location of the 
maximum of the peak. Overall, although the MICE catalogue 
used does not mimic the characteristics of our SDSS samples, 
we can be confident that neither redshift- space distortions nor 
a radial selection function (when it is taken into account in the 
calculation) should bias significantly our results. 

4. SDSS samples used 

We used data from two diff'erent samples of the latest data re- 
lease (DR7) of the spectro scopic SPSS. On o ne side, we used 
the 'Main' galaxy sample (IStrauss et al.ll2002l) as mass tracers 
for reconstructing the overall density field ^(x). On the other, we 
used the LRGs as tracers of the central oyer-densities, and there- 
fore used them as the selected centres x^^ to compute B(R, s). 

Luminous Red Galaxies were selected by the SDSS team 
using several colour and magnitude cuts to obtain a highly bi- 
ased sample reaching high redshift dEisenstein et al .11200 li) . The 
galaxies selected in this way are known to reside ne ar the centres 
of massive dark matter haloes (IZheng et al.ll2009h and are thus 
adequate tracers for the centres of baryon acoustic structures. We 
applied an extra cut in the K-corrected, evolved, g-band absolute 
magnit ude: -23.2 < Mg < -2 L2, as in the previous BAO anal- 
ysis by lEisenstein et al.l (l2005h . This results in an approximately 
volume-limited sample in the redshift range 0.15 < z < 0.30. 

'Main' galaxies in the SDSS constitute a much denser 
sample, and are therefore more suitable to map small density 
changes such as BAO shells. We us ed the 'Main' sam ple from 
the Value-Added Galaxy Catalogue (iBlanton et al.ll2005 ), which 
constitutes a magnitude limited sample in the r band, with r < 
17.6. We applied an extra simple cut, Mf < -20. 

For the conversion of angles and redshifts into co-moving 
distances, we used a fiducial cosmology with the parameters 
Qm = 0.25, Qa = 0.75. In all our analysis we use distances 
in units of Mpc, so that they do not depend on the specific 
value of h. We converted the distribution of the 'Main' galaxies 
into a density field ^(x) binning it into a grid with cubic pixels of 
3 Mpc side. We corrected for the selection eff'ects by weight- 
ing each galaxy by the inverse of the average density at its red- 
shift. As explained below, we performed some tests by slightly 
changing this weighting scheme. Although this weighting may 
not be optimal, it should not aff'ect significantly our results, given 
that the wavelet method does not depend on the local background 
level. We used the density field constructed in this way for the 
calculation of the BAOlet coefficients following equation ([3]). 

In our calculations, we could only use the region in which 
these two samples overlap, which corresponds to the redshift 
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Fig. 4. The BAOlet statistic B for the MICE simulation when 
some observational effects are taken into account. In the top 
panel, we show the B(R, s) obtained when redshift- space dis- 
tortions are included in the simulation. In the bottom panel, we 
show the result when a radial selection function is applied to the 
halo catalogue. In both cases, the contours are drawn at steps of 
5 in B. Solid contours correspond to the results with the obser- 
vational effects included. The dashed contours correspond to the 
original real-space result without selection, i.e., they are identi- 
cal to those in the bottom panel of Fig.[3l 



limits 0.15 < z < 0.26. To minimize border effects in the B(R, s) 
calculation, we defined a buffer region of rbuff = 115 Mpc 
from any of the borders of the 'Main' sample volume. We used 
as centres only the LRGs in the inner volume. This allows us 
to use the density field, as traced by the 'Main' sample galaxies, 
from z > 0.09. In order to minimize angular selection effects and 
border eff'ects, we use a compact area of the sky where the angu- 
lar completeness is nearly uniform. This area c overs 5511 deg^ 
and is defined, in the SDSS survey coordinates (IStoughton et al.l 
2002), by the limits -31.25° < t] < 28.75°, -54.8° <A< 51.8°. 
This results in finally using the density field in a volume of 
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2.2 X 10^ h-^ Mpc^ as traced by A^Main = 198342 galaxies. The 
number of LRGs used as centres (avoiding the buffer region) is 
A^LRG = 1599. 

In Fig. [3 we show a slice of this survey showing both the 
'Main' galaxies and the LRGs. We show how, given the buffer 
used, the LRGs used as centres are located only in an inner vol- 
ume of the larger 'Main' sample. To illustrate the idea of the 
method, we show a zoom around a given LRG galaxy. Even for 
this single centre, a slight over-density of 'Main' galaxies is seen 
at the radii of 105-1 10 /z"^ Mpc. 

As the structures we look for are huge, with radii about 
100 /z"^ Mpc, we have to consider the effect of assumed cos- 
mology (different comoving distances) on our result. In order 
to estimate the distance differences, we compared the distances 
in our adopted MICE cosmology (Qm = 0.25, D a = 0.75) with 
these in the WMAP 7-year cosmological model (iKomatsu et al.l 
120111) . = 0.271, Qa = 0.729. We fixed the redshift dif- 
ference Sz = 0.07 that corresponds approximately to our shell 
diameter of 200 Mpc, and found that this gives distance dif- 
ferences of only a 0.3 and 0.8 per cent at the near and far borders 
of our sample (the MICE distances are larger than the WMAP7 
ones in each case). So, for our nearby volume, the effect is small, 
and does not affect our results given that the statistical uncertain- 
ties are much larger (see next Section). However, this effect will 
be significant for deep samples. 



5. Results for the SDSS samples 

We performed the calculation of B(R, s) for the SDSS in an anal- 
ogous way to the case of the MICE simulation, using the sam- 
ples defined in Section |4l Our results are shown in Fig. [6l As 
above, we mask the region R < 2s. As we are not introducing 
any border correction when calculating the B(R, s) statistic, we 
also mask the region corresponding to the values R > rbuf - s. 
Values obtained at a those large values of R could contain some 
spurious signal, as the calculation of Wr^s would rely on the den- 
sity field in regions outside of the survey boundaries. 

The resulting B(R, s) map is qualitatively very similar to 
that expected, either using an analytical ACDM model (Fig.O, 
or the MICE simulation (Figs. [3] and |4]). This is an indication 
that the observed pattern does not originate from spurious fea- 
tures in the SDSS but is closely related to the large scale struc- 
ture and more specifically the BAO. B(R, s) attains a maximum 
at T^max = 116/z"^Mpc, ^-max = 36/z"^Mpc. This maximum 
is clearly related to the characteristics of the BAO structures 
present in our samples. We studied the robustness of this result 
by changing the weighting scheme applied for the construction 
of the density map (see Section |4]). We did so by capping at dif- 
ferent maximum values the possible weights associated to each 
galaxy, and repeating the calculation of B(R, s) in each case. The 
results were qualitatively similar, obtaining a peak in B(R, s) in 
all cases. However, the position of the peak changed in each case, 
with maximum changes of the order of ±5 Mpc in T^^ax, and 
±10 Mpc in ^max- Therefore, the difference between the po- 
sition of the peak obtained from the SDSS data and that given 
by the MICE simulation is not significant. In any case, we can 
not use the scale and the width of the observed maximum of 
B(R, s) as direct estimates of the radius or width of the shells, 
specially given that our analysis of the possible observational 
biases (Fig.|4]) was only qualitative. 

In order to assess the significance of the BAO detection with 
this method, we focused on the value of B(R, s) obtained at the 
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Fig. 6. The BAOlet statistic B calculated for SDSS data as a func- 
tion of the parameters (R, s). The bottom panel shows the results 
in the full parameter space considered, where we sampled both R 
and s at intervals of 1 /z"^ Mpc. We mask two areas, at the upper 
right and left corners, where our results are not reliable (see de- 
tails in the text). The contours are drawn at steps of 5 for 5 < 
(dotted), B = (solid), and B > (dashed). The top two panels 
show cuts at the arbitrarily chosen values s = 36 Mpc (top) 
and s = 20 Mpc (middle), marked with grey horizontal lines 
in the 2D panel. In these panels, the blue line is B(R, s), while 

the green line and the red band show the mean (B ) and 
interval for the Monte Carlo realizations of random centres. We 
obtain a clear significant peak at different values of s, with a 
maximum for R = 116 Mpc, ^ = 36 Mpc. 

maximum, ^^ax = ^(^max, '^'max) = 22.9 ± 3.'fl. A more thor- 
oughfull analysis would model the B(R, s) statistic in the full 
parameter space. However, given the large co variances between 
measurements at different values of (R, s) we do not expect a 



^ This error in 5max is obtained from the variance of the coefficients at 
the A^LRG different LRGs. However, our significance test is independent 
of this error value. 
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Fig. 5. The SDSS catalogues used. We show a 6 degree thick slice drawn from the catalogues used in our analysis. The red points 
are the LRG galaxies, which form a nearly volume-limited sample at z > 0.15. The 'Main' galaxies of the SDSS are depicted in 
blue. As shown, we use the 'Main' galaxies from z > 0.09, and we use only the LRGs in an inner volume, allowing for a buffer 
region of rbuff = 175 /z"^ Mpc from any of the borders of the 'Main' sample volume. For our analysis, we use the samples covering a 
total area of 551 1 deg^. The radius of a typical BAO shell is shown as a segment. At the top insets we show two orthogonal slices of 
width 20 Mpc, centred on a particular LRG. The BAOlet coefficient map W/?,^(x) has a large value at the position of this centre 
(for 7?max, '5'max), and thus we expect to find a strong BAO signal. The two circles have the radii of 100 and 110 Mpc respectively. 
This is a single BAO shell where the over-density can be appreciated by eye at the right scale. 



large diff'erence from the simple case we consider. We will as- 
sess the probability of finding such a maximum in the case in 
which there are not baryon acoustic structures present in our 
sample. We model this null hypothesis by using randomly dis- 
tributed centres for the calculation of BiR, s) in equation ([5]), 
instead of LRGs. Even using the Wr^s{^) coefficients from the 
observed density field (traced by SDSS 'Main' galaxies), the ex- 
pected value of B{R, s) in this case is (see equation |4l), and 
we expect to obtain a significantly higher signal in the data. In 
this way, we are testing the null hypothesis that, either there are 
not shell-like structures in the density field traced by the 'Main' 
sample, or these shell-like structures are not found preferentially 
around LRG centres. In either case, that would mean that there 
are not B AO-like structures present in our sample. 

To perform the significance test, we generated 10^ random 
realizations of a Poisson process, with the mean number of 
points A^LRG, in the same volume as the LRGs considered in the 
calculation (i.e. taking into account the buff'er zone). For a real- 
ization 7, we use the generated points as our centres ^ to com- 
pute the B{R, s) statistic following equation ([5]), using the Wr^s{^) 
coefficients obtained from the data. We can then obtain the mean 
value B (R, s), and the standard deviation cr^^(R, s) of the 

— MC 

Monte-Carlo realizations of the centres. We show B (R, s) and 
a band of lcr^^(R, s) around it in the top panels of Fig.O 

We now calculate our signal-to-noise ratio at the maximum 
as SNRjnax = Bmax/[cr^^(Rmax,Sjnax)] = 6.60, and assess the 
probability of finding such a large value of SNRj^ax anywhere in 
the parameter space for the Monte Carlo realizations. We used 



SNRjnax instead of directly using ^max because for some regions 
of parameter space, specially at low s, cr^^(R, s) is extremely 
large. Therefore, if we used ^max, we would need to arbitrarily 
restrict the parameter space studied, thus introducing a possible 
a posteriori bias. When using SNRj^ax we sample the full param- 
eter space considered in the calculations (as shown in Fig.O. We 
computed the maximum value of SNR for each realization j in 
the full (R, s) range, SNR^^^^\ The distribution of the values of 
SNR^^^^'' is shown in Fig. [71 where it is compared to the value 
of SNRjnax obtained in the real data. We found that only one of 
the realizations Rave a value of SNR^^^^^ larger than SNR^^^^. 
Thus, the probability of obtaining a maximum with such a large 
SNR in the absence of baryon acoustic structures (our null hy- 
pothesis) is p ^ 10"^, equivalent to a ~ 4 Act detection in the 
Gaussian case. 

However, we should stress here that the significance found in 
this work can not be compared directly to other detection levels 
found in the literature, as it has been stated in the introduction. 
In particular, we are not comparing our results with an analytical 
no-BAO model of B(R, s) (such as that shown in the top panel 
of Fig. O, since to do so would require the detailed modelling of 
all the selection eff'ects aff'ecting the two samples used. 

As explained in Section [2l we can extract more informa- 
tion about the BAO phenomenon in our samples making fur- 
ther use of the BAOlet coefficient maps Wr^s(^). Here, we use 
the coefficient values at the positions of the LRG, for the pa- 
rameters Rmax, Sjnax, which corrcspond to the characteristics of 
the BAO shells present in our samples. In this way, the values 
^max = ^/?max,w ^ measure of how strong is the signal com- 
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SNR^ 



Fig. 7. Histogram showing the distribution of the maximum 
SNR values obtained, in the full (R, s) space, for the 10^ Monte 
Carlo realizations of Poisson-distributed centres (S NR^^^^^). 
This histogram has a mean of 3.72 and a standard deviation of 
0.46. We show as a dashed vertical line the value obtained from 
the data (using the LRGs as centres), SNRj^ax = 6.60. Only one 
of the Monte Carlo realizations give a maximum value larger 
than SNRjns,^. 



ing from a BAO shell around a given point, and in particular, a 
given LRG. Therefore, using Wmax we can localise in configu- 
ration space the regions of the volume covered by our samples 
where the BAO signal is mostly coming from. 

We illustrate this idea in Fig. O where we plot a two- 
dimensional projection of the distribution of the LRGs used as 
centers in our analysis, showing also the value of l^max for each 
of them, following a color gradient. The highest values of Wmax 
correspond to the red points in the plot. 

In Table [H we provide the 10 LRGs used as centers 
with the larger values of Wmax- The whole catalogue of the 
A^LRG = 1599 LRGs used as our centers, and the value of 
^max obtained fo r each of them can be found at the web page 
http : //www . uv . es/martinez This catalogue could be used 
to study the relation of the BAO signal at a given LRG to its 
properties or the environment. It could also be used to make a 
selection of LRG centres with high signal, and use them to re- 
fine the measurements of the BAO characteristics. 

As an illustration of this later use, we show a simple way 
to study the overall properties of the BAO structures, its shape 
and scale. We select those centres which we know that present 
a prominent acoustic feature, i.e., those for which W^ax > 
This leaves us with A^^ = 809 centres. In order to improve the 
signal-to-noise in this illustration for studying the BAO struc- 
tures, we stacked together the 3D density maps around the Nr 
selected LRGs. In doing so, we kept the line-of- sight direction 
aligned for all the centres, as this direction will define the pos- 
sible anisotropics in the distribution. We show a 3D view and 
a 2D cut of this stacked density map in Fig. [9j Thanks to this 
selection the characteristic elements of the BAO are amplified: 
on the one side, a central bump with high density, corresponding 
to the massive halo traced by the LRG, and on the other side, 
the shell surrounding it at a scale of ~ 109 Mpc, showing a 
fainter over-density. We also observe the anisotropic nature of 
these structures. This is a combination of the fact that we have to 



work in redshift space, and of the redshift-dependent selection 
function for 'Main' galaxies. 

A simpler view can be obtained by calculating the average 
radial density profile p(r) around the A^^ centres. The result- 
ing profile, shown in Fig. [TOl has the same features as the 3D 
view: a high bump at short scales, and a clear peak at about the 
acoustic scale, with a maximum at rmax = 109.5 ± 3.9 Mpc. 
The error in r^ax was estimated using bootstrap realizations 
(lLuptonll993D . This scale gives the radius of the baryon acoustic 
shells, and it is therefore a good estimator of the acoustic scale 
in the sample. We also show in Fig. \W[ the radial profiles re- 
stricted to diff'erent regions of the sphere, to better characterize 
the anisotropy of the distribution. We define two cones with a 
width of 45° with respect to the line-of sight in each direction 
(we call these 'near' and 'far' regions), and a 'transverse' region 
covering the belt between the cones. We obtain for each of these 
regions qualitatively similar results. As expected, we see how 
the 'near' and 'far' subsamples are more strongly aff'ected by 
observational eff'ects, such as redshift-space distortions, which 
are more severe along the line of sight. In contrast, the result for 
the 'transverse' subsample matches, within the errors, that for 
the full sphere. It is interesting to note that the value of rmax is 
slightly larger for the 'far' sample than for the 'near' one. 

It is worth to emphasize that this approach would be impos- 
sible with any statistical BAO detection method used this far, 
since the spatial localization of the shells is completely lost in 
the correlation function or in the power spectrum, while the lo- 
cal nature of the wavelet approach has allowed us to identify 
the positions of the most representative structures in our sample. 
Moreover, we are measuring the acoustic scale at positions se- 
lected for their low contamination from other structures, which 
is not the case when averaging over the full sample. In this way, 
we maximize the BAO signal, while minimizing the eff'ect of 
signals coming from diff'erent large-scale structures. 

6. Discussion and conclusions 

In summary, we have designed a new method for the detection of 
baryon acoustic oscillations in the galaxy distribution and for the 
localization, in configuration space, of the structures responsible 
for them. This method is based on the use of a specially designed 
wavelet applied directly on the density field. Our approach also 
relies on the use of two diff'erent tracers: one for the overall den- 
sity distribution, and the other for the central overdensities of the 
baryon acoustic structures. 

After testing the method with simulations, we applied this 
method to the detection of baryon acoustic structures in a sam- 
ple drawn from the SDSS. In this case, we used galaxies from 
the 'Main' catalogue to trace the overall density field, and galax- 
ies from the LRG catalogue to trace the location of massive dark 
matter haloes. We clearly detect BAO in the sample providing a 
confirmation of the detection obtained previously using general 
two point statistics (the power spectrum and correlation func- 
tion). In fact, our approach provides an independent method for 
the detection. Finally, we showed how this method allows us, 
through the use of Wmax(x), to localize in configuration space the 
actual structures responsible for the BAO signal obtained. This 
is a consequence of using a wavelet acting directly on the den- 
sity field. We illustrate the utility of this approach by showing 
the density distribution stacked around a set of centres known to 
show the BAO feature given their T^max value. 

Recent works have p roposed alternative methods to s tudy the 
BAO based on wavelets (IXu et al.l2010l:lTian et al.ll201 ih . In par- 
ticular, iTian et slL use a Mexican hat wavelet function with two 
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Fig. 8. A two-dimensional projection of the distribution in redshift space of the LRGs used in the analysis (i.e. inside our buffer 
region). The data showed corresponds to a slice of a width of ^ 8° in declination, which contains ~ 40% of the LRGs. The colour 
and size of each point corresponds to the value of Wmax at this LRG position. 



SDSS object name 



a (deg) 6 (deg) 



SDSS J141746.20+184733.0 214.44254 18.79250 0.19872 0.00016 517.84 

SDSS J121858.41+380813.6 184.74341 38.13714 0.18974 0.00018 436.68 

SDSS Jl 12430.27+415557.3 171.12613 41.93260 0.19433 0.00019 419.45 

SDSS Jl 12355.53+423816.5 170.98140 42.63793 0.19404 0.00020 414.63 

SDSS Jl 12352.72+424542.4 170.96968 42.76178 0.19469 0.00018 414.63 

SDSS J122935.13+384636.4 187.39640 38.77680 0.18686 0.00016 413.47 

SDSS Jl 12535.99+412608.3 171.39998 41.43564 0.19288 0.00019 401.91 

SDSS J104501. 94+362944.3 161.25810 36.49566 0.15938 0.00015 399.06 

SDSS J140443.3 1+264439.2 211.18047 26.74424 0.15854 0.00016 396.19 

SDSS J142031.28+21 1700.4 215.13036 21.28346 0.19232 0.00020 395.40 

Table 1. Table containing the basic characteristics (J2000 sky coordinates, redshift and redshift uncertainty) of the 10 LRGs with 
the larger values of the BAOlet coefficient at the maximum Wmax- The full table is available at http : //www . uv . es/m artinezj 



parameters, conceptually similar to ours. They use it to search 
for a peak in the two point correlation function of the 'Main' 
SDSS sample, obtaining a detection with a /7-value p = 0.002 
(equivalent to S.lcr in the Gaussian case). As in our case, this 
shows the utility of using the 'Main' sample to reduce the shot 
noise in the calculation and to obtain significant detections. 
However, these works apply the wavelet to the measured two 
point correlation function, instead of directly to the density field. 
In this way, the use the capabilities of the wavelets to character- 
ize accurately the BAO signal (in terms of radius and width), but 
they are not able to get any information about the localization of 
these structures in space. 

The use of wavelets directly on the density field isolates valu- 
able information about the baryon acoustic structures that is hid- 
den in the standard two point statistics. In particular it gives us 
information, through the coefficients Wr^s(^), to localize regions 
in the sampled volume giving the largest or lowest signal. We 
expect that this new method for studying BAO will be of much 
use for ongoing or pla nned surveys, such as the WiggleZ Survey 
(iDrinkwater et al.ll2Q10i) . the Baryon Oscillation Spectroscopic 



Survey (BOSS, 'Eisen stein elaP l201lh . or the Physics of th e 
Accelerating Universe (PAU) Survey (iBemtez et al.l l2009l) . 
which will cover a much larger volume than studied here, and 
will explore higher redshifts. 
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Fig. 9. Stacked 3D and 2D density field. In the top three pan- 
els, we show the density field after stacking the A^^ centres with 
^/^max,w ^ ^- show surfaces encompassing the regions 
above diff'erent thresholds in density after an isotropic Gaussian 
smoothing with cr = 10 /z"^ Mpc. The density threshold de- 
creases from top to bottom, with values of^ = 1.24, 1.18, 1.13 
is the density relative to the average density of the sample). We 
show only the bottom half of the density field for clarity. It can 
be seen that the acoustic shell appears clearly around the cen- 
tral over-density at the detected horizon scale. A 2D slice of this 
density field is shown in the bottom panel. Here, the dotted line 
is a circle whose radius corresponds to the one we measure for 
the BAO shells, r^ax = 109.5 Mpc, and the arrow marks the 
direction of the line of sight. 
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Fig. 10. Average radial profiles. We show the radial density pro- 
file averaged over the Nr centres with Wr^^^^s^^^x ^ ^- P^^t 
S(r) = where po is the average density of the sample. The 
continuous line with red error band shows the radial profile for 
the full sphere. We also show the error bands for the radial pro- 
file restricted to regions of the sphere, as defined in the text: the 
'near' region (green), the 'far' region (blue), and the 'transverse' 
region (orange). The arrow signals the location of the maximum, 
^max = 109.5 /z"^ Mpc. The error band corresponds to the l-cr 
dispersion of 1000 bootstrap realisations. The profiles were es- 
timated using a kernel of width h = 4h~^ Mpc in the radial 
coordinate, but similar results are obtained when using slightly 
diff'erent widths or shapes of the kernel. 
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